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Abstract.  A new  approach  for  calculation  of  resonant  state  parameters  is  developed.  The  method 
proposed  allows  to  solve  different  scattering  problems,  such  as  scattering  and  capture  probability 
as  well  as  calculations  of  shifts  and  widths  of  energy  levels.  It  has  been  applied  to  the  problem  of 
resonant  states  induced  by  impurities  in  the  barrier  of  quantum  wells  and  by  strain  in  uniaxially 
stressed  germanium. 


Introduction 

Quasistationary  (or  resonant)  states  have  been  studied  in  atomic  physics  very  well.  Semi- 
conductors are  other  systems  where  resonant  states  play  a significant  role  in  physical  pro- 
cesses. Such  states  appear,  for  example,  in  gapless  semiconductors  when  doped  by  shallow 
acceptors.  The  system  of  special  interest  is  uniaxially  strained  germanium  where  the  gen- 
eration of  THz  radiation  has  been  achieved  [ , ]. 

Here  we  suggest  a new  method  for  calculating  the  parameters  of  resonant  states  and 
the  probability  of  resonant  scattering,  capture  and  emission  of  carriers.  The  approuch  is 
based  on  the  configuration  interaction  method  which  was  first  introduced  by  Fano  | ] in 
the  problem  of  autoionization  of  He.  The  main  idea  is  to  choose  two  different  hamiltonians 
for  the  initial  approximation  one  for  continuum  states  and  the  other  for  localized  states. 
The  method  is  applied  to  resonant  states  induced  (i)  by  impurities  in  the  barrier  of  quantum 
wells  and  (ii)  by  shallow  acceptors  in  Ge  under  stress. 


1.  Resonant  states  induced  by  localized  states  in  barriers 

We  will  demonstrate  the  general  idea  by  applying  it  to  the  system  consisting  of  a quantum 
well  (QW)  and  one  impurity  in  the  barrier.  The  full  hamiltonian  is  given  by 

- h 2 

H = - — A + V(z)  + Vd(r-r0),  (1) 

2 m 

where  the  potential  of  the  QW  V (z)  is  shown  in  Fig.  1,  IQ  is  the  defect  potential  and 
ro  = (0,  0,  zo)  is  the  position  of  defect. 

As  an  initial  approximation  for  the  wave  function  of  the  localized  state  induced  by  the 
impurity  we  use  the  solution  of  the  equation: 


2m 


A + Vd(r-  r0) 


<P( r - ro)  = Eo(p(r  - r0). 


(2) 
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Fig.  1.  The  parameters  of  the  well  and  impurity. 


The  wave  functions  of  continuum  states  i/'kl  r)  satisfy  the  following  equation: 

lAk(r)  = Ek  tAk(r). 


h2 

- — A + V(z) 
2m 


We  are  considering  a QW  with  one  energy  level  only,  so  that 

Ek  = -AV  + E\  +Sk, 


(3) 

(4) 


where  E\  is  the  space  quantization  level  and  ek  = fi2k2 /2m  is  kinetic  energy  of  the  2D 
motion.  So  we  can  write  for  V'k(i'): 

tAk(r)  = -^<P(z)eikp,  (5) 

where  5 is  a normalizing  square. 

N ow  we  consider  the  problem  of  scattering  of  the  in-plane  moving  carrier  by  the  impurity 
in  the  barrier.  Following  Dirac,  we  construct  the  wave  function  in  terms  of  scattering  theory 
in  the  following  form: 


'l'k(r)  = ^k(r)  +ak<p(r  - r0)  + V — — ^k'(r),  y -*■  0.  (6) 

^ ek  - Sk'  + iy 


As  the  presence  of  one  impurity  does  not  perturb  the  continuum  spectrum  significantly, 
'l'k(r)  should  correspond  to  the  energy  Ek.  Solving  the  Schrodinger  equation  with  the  full 
hamiltonian  (1)  for  'I'ktT)  one  obtains  the  following  expressions  for  a k and  tkk': 


= _1 Vk t 1 Vk  z* 

ak  VsEk-(E0  + AE)  + ir/2’  kk'  S Ek  -(E0  + AE)  + ir/2 

The  energy  shift  A E and  the  width  T /2  of  the  resonant  level  are  given  by: 


4£  = s“a^/‘,2‘'z^'+a^'’/ 

\ f d2k'Z*,VvS(Ek  - Ek>). 

Here  the  matrix  elements 


d2k' 


Z*,v W 
Ek  - Ek' 


(7) 

(8) 

(9) 


Vk  = '/S(<P\VM k),  WU  = Vs(<p |^k>,  5 = (<p\v(z)\<p),  Zk  = VS(<p\ V(zMk),  (10) 


are  introduced. 
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The  probability  of  resonant  elastic  scattering  Wkk'  of  2D  carriers  and  the  capt  ure  prob- 
ability Wkr  are  given  by: 


Wkk'  = — |tkk'|2<S>(Di  —«*')»  Wkr  = |Ok|2-  (11) 

ri 

respectively.  Both  probabilities  contain  the  same  resonant  denominator. 

The  resonant  scattering  should  be  introduced  into  the  kinetic  equation  when  one  solves 
the  problem  of  the  2D  carriers  distribution  function  under  an  electric  field  applied  in  plane 
of  the  quantum  well.  This  scattering  affects  the  distribution  function  /k  of  hot  carriers, 
which  is  connected  with  the  population  of  impurities  in  the  barrier  fr  [ ]: 

fr  = YJ  Wkr/k-  (12) 

k 

2.  Resonant  acceptor  states  in  uniaxially  strained  germanium 


Tetraherally  coordinated  semiconductors  (eg.  GaAs,  Ge,  Si)  have  a fourfold  degenerate 
top  of  the  valence  band.  When  strained,  the  top  of  valence  band  is  split  into  two  doubly 
degenerate  states.  The  ground  state  of  an  acceptor  shows  the  same  behavior  under  uniaxial 
stress.  At  some  critical  value  of  the  stress  — when  the  splitting  is  larger  than  the  acceptor 
binding  energy  — one  of  the  split  levels  is  shifted  into  the  continuous  spectrum  of  the  other 
valence  subband  and  becomes  resonant.  An  effective  optical  transitions  between  resonant 
and  localized  states  of  the  same  impurities  can  take  place.  If  the  electric  field  is  strong 
enough  an  electric  impurity  breackdown  occurs  and  practically  all  localized  impurity  states 
become  depopulated.  Now  capture  and  emission  processes  lead  to  an  effective  population 
of  resonant  states.  This  may  cause  an  intracenter  population  inversion  that  is  the  basis  for 
THz  generation  [ , ]. 

Resonant  acceptor  states  were  considered  in  [ ] by  using  the  Dirac  approach,  which 
requires  choosing  an  initial  approximation  hamiltonian  giving  localized  states  overlapping 
with  the  continuous  spectrum.  The  approach  in  [ ] applies  for  large  stresses  and  for  small 
quasimomenta  but  it  fails  for  the  region  of  the  spectrum  where  resonant  states  are  present. 

Using  our  new  approach  we  will  consider  the  ground  resonant  state  induced  by  shallow 
acceptors  in  uniaxially  strained  p-Ge  along  [00 1 ] direction. 

As  the  initial  approximation  for  localized  states  we  choose  the  diagonal  part  of  the 
Luttinger  hamiltonian  and  the  Coulomb  potential  of  an  acceptor.  For  continuum  states 
we  use  the  eigenfunctions  i/k  (r)  of  the  full  Luttinger  hamiltonian  for  free  holes  in 
cylindrical  approximation.  Following  the  procedure  from  the  first  section  we  are  looking 
for  wave  functions  in  the  form: 


*f/2  = ^ + «k±1/2’+3/V3/2(r)  + af/2-3/ V3/2(r) 


+E 


±1/2, +1/2 

'kk' 

£k  - «k'  + i Y 


<y2  + E 


±1/2,— 1/2 

'kk 

«k  - «k'  + i y 


1/2 


(13) 


We  carried  out  calculations  taking  into  account  non-resonant  scattering  by  the  Coulomb 
potential  in  first  order  only.  The  probability  of  capturing  holes  with  momentum  projection 
+ 1 /2  Wkr  and  the  probability  of  elastic  resonant  scattering  Wkk'  are  defined  now  by 


Wkr  = 


+ 1/2, +3/2 


+ 


+1/2, -3/2 


Wkk'  = 


2 71 


. + 1/2, + 1/2 

kk' 


+ 


+1/2-1/2 

kk' 


<$(£k  ^ £k') 


(14) 


(15) 
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Fig.  2.  The  position  (a)  and  width  (b)  of  resonant  state  as  a function  of  stress  applied  in  [001] 
direction. 


The  expressions  for  the  case  of  momentum  projection  — 1 / 2 are  similar.  The  results  of 
calculations  of  the  resonant  level  shift  and  level  width  are  given  in  Fig.  2. 
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